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Using a simple analytic approach, we study the universal properties of second-order phase 
transition in holographic superconductor models. We explore a general model in arbitrary 
dimensions in which the condensation occurs via the Stiickelberg spontaneous symmetry 
breaking mechanism. All the possible second-order phase transitions and their universal 
characteristics can be identified analytically. The relationship between the critical tem- 
perature and charge density is generic, and the critical exponents can be greater than the 
|' typical mean field value 1/2. In addition, the related numerical factors can also be computed 

D , qualitatively. 

£SJ \ §1. Introduction 

In the past years, the anti-de Sitter /conformal field theory (AdS/CFT) cor- 
respondence^ 1 ^ was intensively applied to study strongly coupled phenomena in 
various physical problems, such as in quantum chromodynamics (QCD), condensed 
matter physics etc. In virtue of the remarkable strong/weak duality, the AdS/CFT 
provides an astonishing advance to deal with certain field theories in the strong cou- 
' pling region via their holographic dual weakly coupled gravitational descriptions. In 

particular, the investigation of the phase transitions in holographic superconductors 
has achieved significant progress, see excellent reviews.®"® The pioneering scheme 
of the gravitational dual description,®'® the so-called abelian Higgs model, consid- 
ered a scalar field and a U(l) gauge field in the planar Schwarzschild-anti-de Sitter 
black hole. A nontrivial profile of the scalar field manifests the condensate of the 
paring mechanism, and the U(l) gauge field exhibits electromagnetic properties in 
superconductors. 

One of the most interesting phenomena in holographic superconductor research 
is the second-order phase transition. The correlation lengths near the critical point 
of second-order phase transitions become divergent, so the systems can be described 
by a scale-invariant theory. In such circumstances, the AdS/CFT is an exceptionally 
appropriate technique to investigate the corresponding universal natures, in particu- 
lar, the critical exponents. It turns out that the critical exponent in the abelian Higgs 
model is 1/2, which is the universal value in the mean field theory. Later, a general- 
ized modef^'E® was proposed in which the condensation occurs via the Stiickelberg 
spontaneous symmetry breaking mechanism. The new features of the Stiickelberg 
holographic superconductor models include the presence of first-order phase transi- 
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tions and second-order phase transitions with nonmean field behavior 

Moreover, some analytic approaches have been proposed to address the univer- 
sal properties of second-order phase transitions in holographic superconductors. One 
suggestion is to consider a simple match of the boundary and horizon solutions of 
both the scalar and U(l) gauge fieldspl see also Refs.^'^ Another analytic ap- 
proach introduces a trial function with a free parameter which that is suggested to 
be fixed by a minimizing procedure,^ see also RefsP-^ In this work, we apply the 
simple matching approach to analytically study the second-order phase transition in 
a general class of the Stiickelberg holographic superconductor in arbitrary dimen- 
sions. We examine all the possibilities that the second-order can occur, and derive 
the corresponding critical exponents. Our results are consistent with the numerical 
analysisP'^ 1 Some related properties in three- and four-dimensional spacetimes 
have been discussed mP^'G^ 

The outline of the paper is as follows. In §2, we review the Stiickelberg holo- 
graphic superconductor models. The analytic study of phase transition is given in 
§3, and explicit models are examined in §4. Some discussions are given in the last 
section. 



§2. Stiickelberg holographic superconductor 

The Stiickelberg holographic superconductor^ 1 '^ i s a simple generalization of 
the abelian Higgs modeP^~® in which the condensation occurs via the Stiickelberg 
spontaneous symmetry breaking mechanism. In the probe limit, the action of the 
model can be decoupled into two parts. The gravitational section provides the back- 
ground of a black hole in the anti-de Sitter (AdS) space providing the temperature 
in the holographic superconductor. The action of this section, in general (d+1)- 
dimensions, is 

S g = J d d+1 xV=g(R-2A), (2-1) 
and the gravitational background is a planar Schwarzschild-anti de Sitter black hole 
ds 2 = -f(r)dt 2 + f~ l {r)dr 2 + r 2 (dx\ + ■■■ + dx 2 ,^) , (2-2) 
where (the radius of AdS is fixed to be unity, i.e., L? = —d(d — 1)/2A = 1) 

f(r) = r 2 (l - ^) . (2-3) 

The parameter rn represents the radius of the black hole, which is determined by 
its mass rfj = M d ~ 2 . The Hawking temperature of the black hole is 

T=±r H . (24) 
In the later analysis, it is more convenient to redefine the radial coordinate as 
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such that the locations of the boundary (asymptotic) and horizon are at z = and 
2 = 1, respectively. 

The second part generically includes a pair of real scalar fields, (ip,<p), and 
one U(l) gauge field, A„, in response to the condensate and conductivity in the 
superconductor. All the fields are assumed as probes and their back reaction to 
gravity can be neglected. The action of this section is 

S m = J d^x^-g (-\f 2 - i(^) 2 - - \rm*P ~ ^) 2 ) • (2-6) 

In this model,® the dynamics is significantly dependent on the choice of the function 
J r (V ; ). The abelian Higgs model indeed corresponds to a particular choice of J 7 ^) = 
t/> 2 . In this paper, we will focus only on the case that J-(ip) is a polynomial of tp. 
The corresponding equations of motion are 

-±=d, (V=gF^) + ?y>W<P ~ A u ) = 0, (2-7) 

=0, (2-8) 

' ^ (v^V) -m 2 ?P-~T'{iP)(d<p-A) 2 = 0. (2-9) 



V~9 

By the gauge freedom, one can assume (p = 0. Moreover, we will only consider 
the scalar potential excitation of the U(l) gauge field, i.e., A = <pdt. Therefore, the 
equations of motion in the gravitational background of (|2-2p reduce to, in terms of 
the coordinate z, 

^"^^-^^^ = 0, (2-10) 

n n d - 1 + z d n , m 2 1 F'(ib)6 2 

d '*-liJ^ d ^-W^* + 2^^ = - <241) 

§3. Analytic study of phase transition 

A simple analytic approach by taking a trivial match of asymptotic and near- 
horizon solutions has been applied to address the properties of the phase transition 
in holographic superconductors.^ It turns out that this simple analysis can really 
capture the universal properties of second-order phase transitions. In this section, we 
are going to reproduce such study for the Stuckelberg holographic superconductor 
model. 

At the boundary, z = 0, the solutions of (f> and tp have the following asymptotic 
forms 

i/) B (z) = D + z x + +D_z x ~, (3-1) 
0B(2) = At--^2^~ 2 , (3-2) 

where 

d± Vd 2 + Am 2 , . 

A± = — n • (3-3) 
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The parameters fi and p are interpreted as the chemical potential and the charge 
density of the dual theory on the boundary. There are two possible condensates of 
the scalar field ip by turning on either the D + or D_ mode. Then, the associated 
vacuum expectation values of the dual operators are 

{0 ± )=r^D ± . (3-4) 

The parameters A± define the dimensions of operators 0±. From the expression 
for A±, the mass of the scalar ip should be m 2 > —d 2 /4, which is consistent with 
the Breitenlohner-Freedman (BF) bound by stability of a scalar field in AdS^+i- A 
typical choice in the literature is m 2 = —2. 

At the horizon, z = 1, the solutions can be expanded as 

<f> H {z) = <M1) " -*) + " z? + • • • > (3-5) 

^ H (z) = ^(1) - ^(1)(1 -z) + \^{z){l -z) 2 + .... (3-6) 

To be regular at the horizon, we should impose an additional condition </>#(l) = 0. 
After that, by solving the equations of motion order by order of 1 — z, the expanding 
coefficients are related as, denoting a = <p' H {\) and b = ^h(I), 

4>^(l) = (d-?,)a-- d aF{b), (3-7) 

9 

Tn 

^(1) = -— 6, (3-8) 

,2,^2 



^(1) = ^ ^ 6 " ^ra^W- (3-9) 



m 2 (m 2 + 2d) . 

"•#/ 

The essential step in any analytic study is to construct a suitable scheme to 
resolve the relations of the coefficients in the boundary and horizon solutions. For 
our simple approach, we just take a trivial manner to match the boundary and 
horizon solutions of 4> and ip, and their first derivatives at a given matching point 
z m in the range of < z m < 1. It turns out that different choices of the matching 
point will give different numerical coefficients, but the universal properties, such as 
the critical exponent, are actually independent of the choice of z m . Let us consider 
a generic choice of z m , and the matching conditions 

4>B(Zm) = (f>H(z m ), (/>' B (z m ) = <t>H{Zm); 1pB{z m ) = ^H(z m ), 1p' B {z m ) = ^' H { z m), 

(3-10) 

give the following relations; here, (D,X) could be either (D + ,A+) or (D_,A_), 

~2 

A* " z ^ 2 = ~~ z ™a ~ ^ [d(d - 3) - F{b)] a, (3-11) 
-(d - 2)^zt 3 = a - ^ [d(d - 3) - F{b)\ a, (3-12) 



2 ~ ~2 

m ° + d + Ad 2 



m 2 (m 2 + 2d)b- -V^(6) 
2r^ 



(3T3) 
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\nJ-l_ 17)2 h z ™ 

XDz m — 



a 2 



m 2 (m 2 + 2d)b - £2-^(6) 



?r z 
zr H 



(3-14) 



where z m = 1 — z m is introduced to simplify the expressions. 

By suitably combining equations (|3-13p and (|3-14p . we can solve D as 

D = x m "^ + 2d 6> (3 .15) 

dzm (Az m + 2z m ) 



and also to get 



where 



A = Aid, m, 3m ) = »V 2 + 2d) + " + *">! . (347) 

^m'^^m "I - 2.ZmJ 

The solution of D shows that the critical exponent of the second-order phase tran- 
sition is indeed determined by the expression of b. Finally, the equation for b can be 
obtained by eliminating a from eqs. (|3-12p and ()3-16|) (taking the minus square root 
to ensure the second-order phase transition) 



1 (,/ - :$)-:„, + ^(b)) = (d ^S 3 ^T- (3-18) 



r 



The solution of b is strongly dependent on the choice of the function J- (if)). In the 
next section, we will explore a general class of interesting models in detail. 

§4. Explicit models 

In this section, we will explore all the possible second-order phase transitions for 
the models that J-{ip) is a polynomial in tp. 

4.1. JXVO = ^™ 

For the simplest case F{ip) = ip n , the master equation ()3-18p for b reduces to 



V d J ^J2A r" H 

From this equation, it is obvious that n = 2 is a very special value in which the 
overall factor ft 1- ™/ 2 disappears. Indeed, this is the only case that the second-order 
phase transition can happen. The value of b should tend to zero near the critical 
point of the second-order phase transition. Therefore, as the temperature approaches 
to the critical value, the overall factor 6 1_n / 2 vanishes for n < 2. The above equation 
then implies p — > 0, which indicates that no condensates occur. For the other case 
n > 2, however, we conclude p — > oo near the critical point. This result seems to 
indicate that a condensate does appear but it is not a second-order phase transition. 
Naturally, one might expect that the phase transition in this case is first order. 
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This expectation is supported by numerical results analyzed in J*** Nevertheless, our 
approach considered in this paper is not capable for studying these two cases. 

The special case n = 2 indeed is just the abelian Higgs model,® '^3 which is well 
studied in the literature, in particular, for dimensions d = 3 and 4. There are nice 
reviews on this model.®"® Generally, the solution for b can be expressed in terms 
of the black hole temperature and critical temperature, T c , as 



d 



[I - (d- 3)z m ] 



^d-1 



fd-1 



1 



T 



d-l 



d-l 



T, 



where the critical temperature is given by 

i (d - 2)4r 3 



d 

Air 



i 

d-l 



(4-2) 



(4-3) 



fA 1 - id - 3)z n 

It is clear that the numerical coefficients in the expressions of b and T c depend on 
the choice of z m . This reveals one weakness of this simple approach: It can only 
give qualitative values for the numerical coefficients. A simple comparison for d = 3 
is given in Table I. However, our approach can capture the universal properties, in 
particular, the relation between the critical temperature and charge density, and the 
critical exponent as 

1/2 



P 



l 

d-l 



D 



1 



T 



d-l 



t; 



d-i 



(4-4) 





Zm = 3/4 


Zm = 1/2 


Zm = 1/4 


Numerical result 


A+ = 2 


T c 
(O) 


0.081 
74.29 


0.104 
71.63 


0.128 
105.28 


T c ~ 0.11 

(O) ~ 144T C 2 (1 - T 2 /T c 2 ) 1/2 


A_ = 1 


T c 
(O) 


0.114 
15.20 


0.142 
11.40 


0.161 
10.05 


T c ~ 0.226^/p 

(O) ~ 9.3T C (1 - T 2 /T 2 ) 1/2 



Table I. "Qualitative" numbers derived by the simple analytic approach for the coefficients in 
the critical temperature and in the condensate of the abelian Higgs model with parameters: 
n = 2,d = 3,m 2 = -2. 



4.2. F(lp) = V 2 + C a l\) a 

In this subsection, we will consider an interesting generalization to the abelian 
Higgs model in which the second-order phase transition can still occur. In this model, 
we consider J-(ip) = i^ 2 + C a ip a . The ip 2 term is needed to ensure the occurrence 
of the second-order phase transition, and the additional term C a ip a with a > 2 is 
a likely generalization. For such a model, the Eq. (|3-18p becomes quite complicated 
and it is impossible to obtain an explicit solution. However, our focus is on the 
properties near the critical point of the second-order phase transition, where the 
parameter b is expected to be small. Thus, Eq. (|3-18|> can be expanded in the order 
of the small parameter b. By keeping up to the leading order of y/b/J-'(b), we have 



C a b a ~' 1 - (d - 3)zm + = " ^ I, - (4-5) 



V2 \ 4^ ; V ' ' d v J V2A 



r d-3 
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It is clear that the dominant term on the left-hand side of the equation is different 
for the two distinct ranges of a: (i) a > 4 and (ii) 2 < a < 4. In addition, the 
particular value a = 4 is a special case. Let us analyze these cases separately. 

For the value of a in the range of a > 4, the b a ~ 2 term is of higher order 
and can be neglected compared with b 2 . Therefore, the critical exponent is exactly 
identical to the abelian Higgs case, namely, the mean field value 1/2. Therefore, the 
additional term C a ip a with a > 4 completely does not change any of consequences 
on the second-order phase transition. 

For the case 2 < a < 4, instead of b 2 , the term b a ~ 2 becomes dominant and 
Eq. (|4-5p reduces to 

This equation allows a second-order phase transition type of solution for b, namely, 

4 T?- 1 ( T d - 1 ' 



,a-2 



ft 



C Q T*- 1 V T r d_1 



with the same critical temperature given in (|4-3p . To ensure that b has a real value, 
some possible constraints on the parameter C a are required; for example, C4 should 
be negative for a = 4. Therefore, the critical exponent in 

d-i 1 

rpa-2 / T^rf-1\^2 

is l/(a— 2), which is always bigger than the mean field result. Note that the exponent 
is discontinuous at a = 2 of the mean field result. 

Finally, a = 4 is a special value and in such a case, we should solve 



1 

V2 



L - (</ - ■■>)-:,„ + ( *p - [1 - (d- 3)z m ]Q ) 6 2 (</ ~ ~ >:r> ^ 



2Z\ rl" 1 ' 



(4-9) 



The coefficient of b 2 should be positive, which introduces a constraint to the param- 
eter C4 as 

< , ri /; m T . (4-io) 

d[l - (d - 3)z m ] 

However, the explicit value of this constraint is not exact. It depends on the choice 
of z m , so our analysis can only give a qualitative result. Obviously, the critical 
exponent is identical to the mean field result, but the numerical factor is modified 



d[l-(d-3)zm\C4T*-i V Tt l 



which again cannot be exactly obtained by our approach. 
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4.3. H^)=^ 2 + EC a ^ ai 

The most general model for a polynomial J-(ip) 1S •^(' l P) = ^ 2 + C ai ip ai with 
aj > cti > 2 when j > i. Again, the term if 2 is necessary for the existence of the 
second-order phase transition. It is not difficult to recognize a general feature: the 
much higher order extended terms are indeed more insignificant at the critical point. 
Therefore, for the study of universal properties, in addition to ip 2 , one only needs to 
keep the lowest order correction term, assuming C ai vp ai . Hence, a complete analysis 
was already given in the previous subsection. 

§5. Conclusions 

In this paper, we analytically study the universal properties of the second-order 
phase transition for a general class of holographic superconductors, the Stiickelberg 
model with a polynomial J-(ip). The desired characteristics can be revealed by a 
simple approach, namely, by matching the asymptotic and horizon solutions of two 
essential scalar fields at an arbitrary point in between. The validity of this approach 
can be understood owing to the presence of a scale invariance at the critical point 
of the second-order phase transition. Firstly, we can classify all the possible cases 
admitting the second-order phase transition. Moreover, we can derive the explicit 
expression of the critical exponent. For most cases, the second-order phase transition 
has mean field critical exponent 1/2. The essential exceptional example is the case 
= tp 2 + C a ip a with 2 < a < 4. In this case, the critical exponent is (a — 2) _1 , 
which is always greater than the mean field value. All our results are consistent with 
the numerical analysis P'EB 

A significant weakness of our simple approach is that it cannot give very exact 
numerical coefficients since those values are dependent on the choice of the match 
point. Another type of analytic approach has been proposed, requiring an eigenvalue 
minimization while introducing a trial function^} It has been checked that this 
approach can give very accurate results for the numerical coefficient in many models, 
including an external magnetic fielcP^J and the Gauss-Bonnet gravityP^J'^ One 
should expect to improve the precision of the numerical coefficients in our results by 
considering other analytic approaches. 
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